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Abstract. It has been demonstrated, using variational methods, that quantum 
vacuum energy gravitates according to the equivalence principle, at least for the finite 
Casimir energies associated with perfectly conducting parallel plates. This conclusion 
holds independently of the orientation of the plates. We review these arguments and 
add further support to this conclusion by considering parallel semitransparent plates, 
that is, (5-function potentials, acting on a massless scalar field, in a spacetime defined 
by Rindler coordinates. We calculate the force on systems consisting of one or two 
such plates undergoing acceleration perpendicular to the plates. In the limit of small 
acceleration we recover (via the equivalence principle) the situation of weak gravity, and 
find that the gravitational force on the system is just A/g, where g is the gravitational 
acceleration and M is the total mass of the system, consisting of the mass of the plates 
renormalized by the Casimir energy of each plate separately, plus the energy of the 
Casimir interaction between the plates. This reproduces the previous result in the 
limit as the coupling to the (5-function potential approaches infinity. Extension of this 
latter work to arbitrary orientation of the plates, and to general compact quantum 
vacuum energy configurations, is under development. 



PACS numbers: 03.70.+k, 04.20.Cv, 04.25.Nx, 03.30.+p 



1. Introduction 

The subject of Quantum Vacuum Energy (the Casimir effect) dates from the same year 
as the discovery of renormalized quantum electrodynamics, 1948, and suggests that the 
assertion that zero-point energy is not observable is invalid. (For a contrary viewpoint 
see Ref. [1].) On the other hand, because of the severe divergence structure of the 
theory, controversy has surrounded it from the beginning. Sharp boundaries give rise to 
divergences in the local energy density near the surface, which may make it impossible 
to extract meaningful self-energies of single objects, such as the perfectly conducting 
sphere considered by Boyer [2]. These objections have recently been most forcefully 
presented by Graham, et al. [3] and Barton [1], but they date back to Deutsch and 
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Candelas OIB]. In fact, it now appears that these surface divergences can be dealt with 
successfully in a process of renormalization, and that finite self-energies in the sense of 
Boyer, may be extracted [ZIIH]. 

But the most troubling aspect of local energy divergences is in the coupling to 
gravity. The source of gravity is the local energy-momentum tensor, and such surface 
divergences promise serious difficulties. As a prolegomenon to studying such questions, 
we here address in section [2] a simpler question: How does the completely finite Casimir 
energy of a pair of parallel conducting plates respond to gravity? (We'll address 
divergences in section [31 ) The question, and its answer, turn out to be surprisingly 
less straightforward than the reader might suspect! (For a complementary view on the 
gravitational effects of Casimir energy, see the contribution to these Proceedings by S 
A Fulling et al.) 



2. Variational method 



2.1. Casimir stress tensor for parallel plates 

Brown and Maclay [9j showed that, for parallel perfectly conducting plates separated by 
a distance a in the ^-direction, the electromagnetic stress tensor acquires the vacuum 
expectation value between the plates 

{Tn = 7diag(l, -1, -1, 3), S, = -T^.nc. (1) 

Outside the plates the value of {T'^'^) = 0. Because there are some subtleties here, let us 
review the argument for the case of a conformally coupled scalar (the electromagnetic 
case differs by a factor of two). Actually, the result between the plates, < 2; < a is 
given in great detail in Ref. [lOj (7 is the conformal parameter): 

= (% + M)diag(l, -1, -1, 3) + (1 - 67)(7(z)diag(l, -1,-1,0), (2) 

where 

1 f°° tt'^ 

Uq = / duK^, u = 7. (3) 

" Utt^ Jo ' 1440a4 ^ ' 

Note that uo is a divergent constant, independent of a, and is present (as we shall see) 

both inside and outside the plates, so it does not contribute to any observable force or 

energy (the force on the plates is given by the discontinuity of {T^z)), and so may be 

simply disregarded (as long as we are not concerned with dark energy). But see below! 

Similarly, the term involving the Hurwitz zeta function, 

9iz) = -Y^[C(4, z/a) + C(4, 1 - z/a)l (4) 
which exhibits the universal surface divergence near the plates, 

is also unobservable (if we disregard gravity) because it does not contribute to the force 
on the plates, nor does it contribute to the total energy, since the integral over g{z) 
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between the plates is independent of the plate separation. Of course, the best way to 
eliminate that term is to choose the conformal value 7 = 1/6. 

Since the exterior calculation does not appear to be referred to in Ref. \10\, let us 
sketch the calculation here: Consider parallel Dirichlet plates at 2; = and z = a. The 
reduced Green's function satisfies 

/ d2 ^ 



dz 



+ ^'\g{z,z') = 6{z-z'), (6) 



where = k"^ — uj"^ = k"^ + The solution for z,z' < is 

g{z, z ) = --e'"^< sinh kz>. (7) 

K 

It is very straightforward to calculate the one-loop expectation value of the stress tensor 
from 

i(T'^'^) = (d^'d"' - Ig'^'d^d'^ Gix, x') - 7(9'^9" - g>"'d^)Gix, x). (8) 
V 2 / x'=x 

After integrating over uj = i( and k, we find the result {z < 0) 

{Tn = «odiag(l, -1, -1, 3) - ^^^:^diag(l, -1, -1, 0). (9) 

This is exactly as expected. The uq term is the same as inside the box, so is just the 
vacuum value, and the second term is the universal surface divergence (independent of 
plate separation), which can be eliminated by setting 7 = 1/6. 

Thus, we conclude that the physical stress tensor VEV is just that found by Brown 
and Maclay: 

{T^"") = u diag(l, -1, -1, ?,)e{z)e{a - z). (10) 
in terms of the usual step function. 



2.2. Variational principle 

Now we address the question of the gravitational interaction of this Casimir 
apparatus [H]. It seems this question can be most simply answered through use of 
the gravitational definition of the energy-momentum tensor, 

= -i I {dx)^5g^''T,, = {dx)^5g,,T^r (11) 

For a weak field, gfj,u = rj^i, + 2h^y (Schwinger's definition [12] of h^^). So if we think 
of turning on the gravitational field as the perturbation, we can ignore ^/—g. The 
gravitational energy, for a static situation, is therefore given by {5W = — J dt 6E) 

Eg = - j {d^)K,T^r (12) 

We can use the gravity- free electromagnetic Casimir stress tensor ( fTOj) . with u now 
replaced by Ed a for the electromagnetic situation. 
We now use the metric [131 E] 

5(00 = - (1 + 25(z) , gij = 6ij. (13) 
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Figure 1. Relation between two Cartesian coordinate frames: One attached to the 
earth {x,y,z), where —z is the direction of gravity, and one attached to the parallel- 
plate Casimir apparatus (C, ?7, x)j where C, is in the direction normal to the plates. The 
parallel plates are indicated by the heavy lines parallel to the 77 axis. The x = x axis 
is perpendicular to the page. 



This is appropriate for a constant gravitational field. (But see below.) Let us consider 
a Casimir apparatus of parallel plates separated by a distance a, with transverse 
dimensions L ^ a. Let the apparatus be oriented at an angle a with respect to the 
direction of gravity, as shown in figure [H Let us take the Cartesian coordinate system 
attached to the earth to be {x,y,z), where, as noted above, z is the direction of — g. 
Let the Cartesian coordinates associated with the Casimir apparatus be {(, 77, x), where 
( is normal to the plates, and 77 and x ^^re parallel to the plates. The relation between 
the two sets of coordinates is 

z = ( cosa + rjsma, y = rjcosa — (sina, x = x- (14) 

Let the center of the apparatus be located at (Co, = 0, x = 0). 
Now we calculate the gravitational energy 

f nn /"^/^ /•Co+a/2 

Eg = {dx^gzT = — gL / drj d( {( cos a + rj sin a) 

J a J-L/2 

= ^^L^ cosaa^o + -^5 (15) 
a 

where i^' is a constant, independent of (q. Thus, the gravitational force per area on the 
apparatus is independent of orientation 

^- = = -TT^c = -^^c, 2;o = Cocosa, (16) 
A Aozq 2a 

a small upward push. Here e = 2ga is a measure of the gravitational force relative to 
the Casimir force. Note that on the earth's surface, the dimensionless number e is very 
small. For a plate separation of 1/im, 

e = ^ = 2.2 X 10-^\ (17) 

so the considerations here would appear to be only of theoretical interest. The effect is 
far smaller than the Casimir forces between the plates. 
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It is somewhat simpler to use the energy formula to calculate the force by 
considering the variation in the gravitational energy directly, as we can illustrate by 
considering a mass point at the origin: 

T^"" = m5(r)5'^°5'^°. (18) 

If we displace the particle rigidly upward by an amount Szq, the change in the metric is 
^hoo = —gSzQ. This implies a change in the energy, exactly as expected: 

6Eg = —mc^ {—g6zo) = mg5zQ. (19) 

Now we repeat this calculation for the Casimir apparatus. The gravitational force 
per area on the rigid apparatus is ^ = — = —g^o again the same result found in 
f|T6|l . which agrees with the second result found by Calloni et al. [13] but is 1/4 that 
found by Bimonte et al. [H], who reproduce the first result of Ref. [13]. Our result is 
consistent with the principle of equivalence, and with one result of Jaekel and Reynaud 

2.3. Alternative calculation 

As in electrodynamics, we should be able to proceed, starting from the definition of the 
field 



6W = J {dx)6T^"'h^^. (20) 

Again, check this for the force on a mass point, with stress tensor given by (fTSl) . so 

6Eg = - J (dr)m {-6r ■ V) (5(r)/i°° = -mSr ■ Vh^°. (21) 

Since = —gz, we conclude — ^ = F = —mgz. 

For the constant field the force on a Casimir apparatus is obtained from the change 
in the energy density 

= -eia/2 - C + Co)^(C - Co + a/2), (22) 
a 

that is, recalling that zq = Co cos a, 

ST'o ^ ^szo^ m - Co - a/2) - 6{( - Co + a/2)] , (23) 
a cos a 

which yields a result identical to (fT6l) [h^^ = —g{(cosa + r^sina)] 

()E F £ 1 C=Co+a/2 

C=fo-a/2 



Adzo A a cos a 
2.4. Metric near the surface of the earth 

However, the above metric (fT3|) . while sufficing for massive Newtonian objects, might 
seem inappropriate for photons. Rather, shouldn't we use the perturbation of the 
Schwarzschild metric, which for weak fields {GM/r <^ 1) is in isotropic coordinates 

m- 

, 2 2GM\ 2,2 2GM\ , 2o 

= - U ) + ( + ^ ^ 
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If we expand this a short distance z above the earth's surface, of radius R, we find 

2GM 



9oo 



R 



+ 2gz 



9ij 



Si 



(26) 



Now, for our Casimir apparatus shown in figure [H each component of the Casimir 
stress tensor contributes with equal weight: 



6E^ 
ASzq 



11 



i22 



^33 



(27) 



T'^x = 0, which is twice the previous result. Note that again the result is 



since T 

independent of a. If instead, we use the second method we have 



1 



STf"" = -fc — (1,-1, -1,3)- 

a cos a 

so, again we get the same result: 



MC-Co + 7T -MC-Co-7T 



6E^ 

5zq 



a 



(dr) 



cos a 



MC-Co- 



(-2^7^) 



f28) 



-A2g£,, (29) 



where z = C, cos a + ?7 sin a. 

We might think we would be able to obtain the same result using the original 
Schwarzschild coordinates, where /iqo = ~gz, hpp = —gz, and all other components of 
hfj^i, are zero. However, now if we use the first method above, the result is proportional 
to 7^00-)-T^P = ^4cos^a, which implies (fortuitously) Bimonte et al.'s earlier result [1 
for a = 0: = —4:g£c cos^ a. 



2.5. Gauge noninvariance 

The reason we get different answers in different coordinate systems is that our starting 
point is not gauge invariant. Under a coordinate redefinition, which for weak fields is a 
gauge transformation of h^^^, [12], h^^, h^^y + d^C,^ + d^^p, where is a vector field, the 
interaction W is invariant only if the stress tensor is conserved, d^^T^^ = 0. Otherwise, 
there is a change in the action, AW = —2 j'{dx)C,udfj,T'^'^ . 

Now in our case (where we make the finite size of the plate explicit, but ignore edge 
effects because L 3> a) 

= ^diag(l, -1, -1, 3)^ (^C - Co + ^ - C + Co) 

X e{r] + L/2)9{L/2 - r])9{x + L/2)9{L/2 - x). (30) 

Thus the nonzero components of d^T^'^ are 

'if 

q^TK = _£ ^s{C - Co + a/2) - 6{C - Co - a/2)] 9..., (31a) 
d.T^"^ = - ^ [6{r] + L/2) - 6{r] - L/2)] 9..., (316) 
S^Tf^x = _ ^ [Six + L/2) - 6ix - L/2)] 9..., (31c) 
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where 9 . . . refer to the remaining step functions. Therefore, the change in the energy 
obtained from AW is 

AEg = ^J dr]dx MCo - a/2,r],x) - ^dCo + a/2,r],x)] 
2F r 

--^y dCdx [UC,-L/2,x)-U(^L/2,x)] 
2£ r 

_ _£ y dC dry [^^(C, V, -L/2) - idC. V, L/2)] . (32) 

2.6. Fermi coordinates 

Since we have demonstrated that the gravitational force on a Casimir apparatus is 
not a gauge- invariant concept, we must ask if there is any way to extract a physically 
meaningful result. There seem to be two possible ways to proceed. Either we add 
another interaction, say a fluid exerting a pressure on the plates, resulting in a total stress 
tensor that is conserved, or we find a physical basis for believing that one coordinate 
system is more realistic than another. The former procedure is undoubtedly more 
physical, but will yield model dependent results. The latter apparently has a natural 
solution. 

A Fermi coordinate system is the general relativistic generalization of an inertial 
coordinate frame. Such a system has been given by Marzlin [l7j for a resting observer 
in the field of a static mass distribution. It is actually a priori obvious that in such 
a system gij is quadratic in the distance from the observer. Thus the "constant field 
metric" is simply the Fermi coordinate metric for a gravitating body, 

ds^ = -(l + 2c/2)dt2 + dr'^ (33) 

Thus, coordinate lengths don't depend on z. The metric ( IT3l) is indeed appropriate, and 
the corresponding gravitational force is therefore given by the result found in that case, 
F/A = -gSc, as in (fT6l). 



2.7. Gauge transformation 

Now we can use the method described in ( l32l) to transform the energy in isotropic 
coordinates to that in Fermi coordinates. We compute the additional gravitational 
energy, in terms of the gauge field which carries us from isotropic coordinates to 
Fermi coordinates, /i^^ = /i^^, + 9^^,^ + d^C,^. Here from and (1551) 

^00 = -9z, hfj = -gz6ij, h^^ = -gz, hfj = 0. (34) 

The gauge field turns out to be 

= Qc^cosa + C??sina^ + fiv,x), (35a) 
^r, = ^9 (Cvcosa + ^r/^sina^ + 9{C,x), (35 &) 

= 7;9{Ccosa + r]sma)x + h{C,r]), (35c) 
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rL/2 rL/2 I 2Sc f<0+a/2 rL/2 I 

/ dr] dx-gcosa{-2Coa) / dC / dx-g cosa{-L)C 

J-L/2 J~L/2 4 a JCo-a/2 J-L/2 2 

Ko+a/2 rL/2 I 

/ d( drj - g {(^ cos a + rj sin a) {—L) 

JCo-a/2 J-L/2 2 



where the functions f , g^ and h are irrelevant. Substituting this into the expression for 
AEg, ([32]), we obtain 

a J-L/2 ' J~L/2 "4 . - , ^ J(o-a./2 

_ 2S, r(0+a/2 .L/2 

a JCo-a/2 J-L/2 

= - AgScCoCosa = -AgS^zo, (36) 

which when differentiated with respect to Zq gives an additional force, — = ^ = 
gSc- When this is added to the isotropic force (127|) . we obtain the Fermi force, 
pi + AF 

= -2gSc + gSc = -gSc = (37) 

as given in (fT6l) . This answer is the second one given in Calloni et al. |T3], but is not 
referred to in the 2006 Bimonte et al. paper [M]. Those authors have now modified their 
analysis to agree with ours |18j . 

3. Rindler coordinates 

We now turn to the consideration of the Casimir apparatus undergoing uniform 
acceleration |19j . Relativistically, uniform acceleration is described by hyperbolic motion 

t = C, sinh r, z = C, cosh r, (38) 
where is the proper acceleration, which corresponds to the metric 

dt^ - dz^ = e'dr^ - de'. (39) 
The d'Alembertian operator takes on cylindrical form 

3.1. Single accelerated plate 

For a single semitransparent plate at ^i, the Green's function can be written as 



9 = ^5(^-0, (42) 



where the reduced Green's function satisfies {k = |k_L|) 

which we recognize as just the semitransparent cylinder problem with m ( = —iu 
and K — > fc. Thus, the Green's function for a single plate is 

1 + /i^i/c(A;^i)i^c(^6) 

where the strong coupling limit, /i oo, corresponds to Dirichlet boundary conditions. 
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3.2. Minkowski- space limit 

If we use the uniform asymptotic expansion (UAE), based on the hmit 

^ ^ oo, ,^1 ^ oo, ^ — ^1 finite , ( = (^i ^ oo, ( finite , (44) 
we recover the Green's function for a single plate in Minkowski space, 

— ■ (45) 



where k = a/P~+~C^, uj = iC- 



3.3. Energy-momentum tensor 



The canonical energy- momentum for a scalar field is given by T^i, = d^cpdi^cp -\- g^i^-^^C, 
where the Lagrange density includes the 5-function potential. Using the equations of 
motion the energy density is 

The force density is given by 

fx = --L9,(v^T^) + Ir^^^Oxg,., (47) 

or 

f^ = -^d^UT^^)~^T^'. (48) 

When we integrate over all space to get the ("coordinate") force (per area), the first 
term is a surface term which does not contribute: 

e 

which when multiplied by the gravitational acceleration g is the gravitational force/area 
on the Casimir energy. Using the expression fH6l) for the energy density, and rescaling 
C = C^) we see that the gravitational force is merely 



/de«/« = - / 7fr„o, (49) 



This result is an immediate consequence of the general formula 

^c = -^/(dr)/^2(^^g(r,r), (51) 
in terms of the frequency transform of the Green's function, 

G{x,x')= r ^e-(*-*')^(r,r'). (52) 
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3.4- Force on single plate 

Alternatively, we can start from the following formula for the force density for a single 
semitransparent plate, 

/« = ^0'95;"^X?-?i), (53) 
or, in terms of the Green's function, 

(54) 

For example, the force on a single plate is given by 



Expanding this about some arbitrary point ^o, with C, = CCo, and using the UAE, we 
get (a is an arbitrary scale to make y dimensionless) 



^ ^ \ I 7^ ^"[^ + t^iiimi)Kdki.)l (55) 



1 Z-"" dyy^ 

which is just the negative of the (divergent) quantum vacuum energy of a single plate. 

3. 5. Two accelerated plates 

For two plates at ^i, ^2, for ^, ^' < ^1, 

g{^ e') = I<K> - /^l^l^l + - /^l/^266^1^2(i^2/l - Kll2) jj^ ^g^^ 

where 

A = (1 + //ieiXl/l)(l + lJi2^2K2h) - Hlli2^1^2llKl, (58) 

and where we have used the abbreviations /„ = hi^^a), I — -^c(^0> ^1 — -^c(^C')) ^^c. 
For e,e'>6, 

(C c'\ T + 1^2^211 + lilll2^1^2hl2{l2Ki - I1K2) , . 

g{t,,t, ) = I<Ky — A A,, (59) 

and for ^i, ^2, for ^1 < < 6, 

_ !^MilJ^2i^KK, + t^lB^Mll^iK, + KI,). (60) 
In the ^0 ^ 00 hmit, the UAE gives, for ^1 < C < 6 (0 = 6" Ci) 

i^gii, O ^ + ^ cosh K{i - a 

_ /fi 
2k 

with 



2kJ 2k V 2kJ 



(61) 



which is exactly the expected result. The same holds in the other two regions. 
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3.6. Force on two-plate system 

In general, we have two alternative forms for the force on the two-plate system: 

^=-(%+%)i/^lnA, (63) 



which is equivalent to 



dCd^k 



From either of these two methods, we find the gravitational force on the Casimir energy 
to be in the ^ oo limit 

T= ^/ dKK^lnAo, Ao = e-2'^"A. (65) 

47r^ JO 

Explicitly, 

1 /-OO 1 H i \ i 

J^=^— dyy' 



967r2a3 Jo ( _M_ > i] ( > A (.y _ i 

\tiia I \iJ,2a I 



1 



1 1 

+ 



Mia /i2(i 



-(4 + ^di+^d2), (66) 



which is just the negative of the Casimir energy of the two semitransparent plates. The 
divergent terms are just the sum of the Casimir energies of each plate separately, (156|1 . 
which serve to simply renormalize the mass/area of each plate: 

^totai = mi + ma + E^x + + = Mi + + 4, (67) 

and thus the gravitational force on the entire apparatus obeys the equivalence principle 

gT =-g{M^^M^ + E,). (68) 

Saharian et al. |20j earlier reached a similar conclusion, but only for the finite part of 
the energy. 

4. Conclusions 

• We have found, after a certain confusion, an extremely simple answer to how 
Casimir energy gravitates: just like any other form of energy, 

^ = -gEc (69) 

This result is independent of the orientation of the Casimir apparatus relative to 
the gravitational field. This refutes the claim sometimes attributed to Feynman 
that virtual photons do not gravitate. 

• Although gravitational energies have a certain ill-defined character, being gauge- 
or coordinate- variant, this result is obtained for a Fermi observer, the relativistic 
generalization of an inertial observer. 
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• This conclusion is supported by an explicit calculation in Rindler coordinates, 
describing a uniformly accelerated observer. This demonstrates, quite generally, 
that the total Casimir energy, including the divergent parts, which renormalize the 
masses of the plates, possesses the gravitational mass demanded by the equivalence 
principle. 

• New developments of this work are in progress, and will be described in part in the 
contribution to this proceedings by K V Shajesh. 
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